Homoclinic Chaos in the Dynamics of a General Bianchi IX Model 
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The dynamics of a general Bianchi IX model with three scale factors is examined. The matter 
content of the model is assumed to be comoving dust plus a positive cosmological constant. The 
model presents a critical point of saddle-center-center type in the finite region of phase space. 
This critical point engenders in the phase space dynamics the topology of stable and unstable four 
dimensional tubes Rx S 3 , where R is a saddle direction and S 3 is the manifold of unstable periodic 
orbits in the center-center sector. A general characteristic of the dynamical flow is an oscillatory 
mode about orbits of an invariant plane of the dynamics which contains the critical point and a 
Friedmann-Robertson- Walker (FRW) singularity. We show that a pair of tubes (one stable, one 
unstable) emerging from the neighborhood of the critical point towards the FRW singularity have 
homoclinic transversal crossings. The homoclinic intersection manifold has topology R x S 2 and is 
constituted of homoclinic orbits which are bi-asymptotic to the S 3 center-center manifold. This is 
an invariant signature of chaos in the model, and produces chaotic sets in phase space. The model 
also presents an asymptotic DeSitter attractor at infinity and initial conditions sets are shown to 
have fractal basin boundaries connected to the escape into the DeSitter configuration (escape into 
inflation), characterizing the critical point as a chaotic scatterer. 



The longtime debate on the chaotic dynamics of 
general Bianchi IX models started with the work 
of Belinskii, Khalatnikov and Lifshitz (BKL) on 
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the oscillatory behaviour of such models in their 
approach to the singularity Q. They showed that 
the approach to the singularity^ — > 0) of a gen- 
eral Bianchi IX cosmological solution is an oscil- 
latory mode, consisting of an infinite sequence of 
periods (called Kasner eras) during which two of 
the scale functions oscillate and the third one de- 
creases monotonically; on passing from one era to 
another the monotonic behaviour is transfered to 
another of the three scale functions. The length of 
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each era is determined by a sequence of numbers 
X s ,0 < X s < 1, s — integer, each of which arises 
from the preceding one by the transformation X s+ i 
— fractional part of 1/X S . From the properties 
of this map it is obtained that the behaviour of 
the model becomes chaotic on approaching the 
singularity (t = 0) for arbitrary initial conditions 
given at to > 0. With the advent of powerful nu- 
merical resources the interest in the behaviour of 
these models revived but - as in the BKL work - 
the procedure has been basically to obtain maps 
which approximate the dynamics of the model de- 
scribed by Einstein's equations and which exhibit 
strong stochastic properties^]. How well these 
discrete maps represent the full nonlinear dynam- 
ics has been subject of much research, particu- 
larly by Berger[|| ^] and RughQ. The interest 
in the chaoticity of Bianchi IX models has been 
mainly focused on the Mixmaster case (vacuum 
Bianchi IX model with three scale factors^), but 
the chaotic dynamics of other Bianchi IX model 
universes has also been discussed in the literature 
(cf. |t]] and references therein, and ||); homo- 
clinic chaos in axisymmetric Bianchi IX universes 
with matter and cosmological constant has been 
treated in JsJ . The chaoticity in the Mixmaster 
dynamics has been object of much dispute in the 
literature (cf. the contributions to the Section 
Bianchi IX (Mixmaster) dynamics in |]l0| , partic- 
ularly Q). Latin et al.[[fl| and Contopoulos et 
al.|12[| have shown the nonintegrability of the Mix- 
master model in the Painleve sense, although the 
question of the generic behaviour (chaotic or not) 
remained unsettled mainly due to the absence of an 
invariant (or topological) characterization of chaos 
in the model (standard chaotic indicators as Lia- 
punov exponents beingcoordinate dependent and 
therefore questionable jl3| Jl(J ) . More recently Cor- 
nish and Levin |l4| proposed to quantify chaos in 
the Mixmaster universe by calculating the dimen- 
sions of fractal basin boundaries in initial condi- 
tions sets for the full dynamics, these boundaries 
being defined by the code associated with one of 
the three outcomes on which one of the three axes 
is collapsing most quickly, as established numeri- 
cally. Their result received a recent critical review 
that nevertheless endorses it|15(]. In the present 
paper we shall give an invariant characterization 
of chaos for the Mixmaster model with a positive 
cosmological constant, but this criterion does not 
work for the case of zero cosmological constant. 

Our purpose in this paper is to examine the dy- 
namics of a general (three scale factors) Bianchi 
IX cosmological model with dust and cosmological 
constant and establish the existence of chaos using 
the criterion of the homoclinic transversal crossing 
of topological tubes that organizes the dynamical 
flow in the model. We show that the dynamics 



of the model is highly complex and chaotic, and 
that chaos has a definite homoclinic origin. The 
phase space of the model is noncompact and the 
presence of the cosmological constant determines 
two crucial facts in phase space: first, the exis- 
tence of a critical point of the type saddle-center- 
center; second, two critical points at infinity corre- 
sponding to the DeSitter configuration, one acting 
as an "attractor" to the dynamics. With respect 
to the latter point, our system has the character- 
istic of a chaotic scattering system with two ab- 
solute outcomes consisting of (i) escape into in- 
flation (the DeSitter attractor) or (ii) recollapse 
to the singularity. The presence of the critical 
point of saddle-center-center type is responsible for 
a wealthy and complex dynamics, engendering in 
phase space topological structures which are ho- 
moclinic to a manifold of periodic orbits, having 
the topology of the three sphere S 3 . The dynam- 
ics has a 2-dim invariant plane allowing to expand 
the dynamical equations about it, in the same pro- 
cedure used to examine the dynamics about a peri- 
odic orbit of the system. The mixmaster universe 
with a cosmological constant will correspond to the 
limiting case of dust equal to zero but, as we will 
discuss, the skeleton of its dynamics is framed in 
the dynamics of the general case. 

The line element of the model has the form 



ds 2 = dt 2 - 



q3\2 



(1) 



where t is the cosmological time and 9 1 = Mfyw 1 , 
9 2 = N(t)uj 2 and 6> 3 = R(t)uj 3 . Here the oj 1 are 
Bianchi IX 1-forms satisfying dui 1 = e lJ a; J A ui k . 
The matter content of the models is assumed to 
be a pressureless perfect fluid, namely dust, with 
energy density p, as described by the comoving ob- 
servers with four velocity orthogonal to the homo- 
geneity surfaces of ((!]), plus a positive cosmologi- 
cal constant A. The dynamics of the three scale 
factors M(t), N(t) and R(t) will be given by Ein- 
stein's equations, which are equivalent to Hamil- 
ton's equations generated by the Hamiltonian con- 
straint 



H 



If M 2 N 2 
~NR Pm ~ MR Pn 

J^Pr + ^PMPN + ^PNPR + ^PMPR 

^L^(i? 4 + Af 4 + N A — (R 2 N 2 ) 2 - (R 2 - M 2 ) 2 

r2 Ar2\2\ 



- (M - N 2 Y) - 2AMNR - 2E Q = 0. 

Here pm , P n and pr are the momenta canonically 
conjugate to M, N and R, respectively. Eq is a 
constant proportional to the total energy of the 
model, arising from the first integral of the Bianchi 
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identities, pMNR = Eq. Hamilton's equations 
have the form 



M = 



dH 

dp 



Pm = 



A I 



dH 
~dM 
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8M 2 N 
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(3) 
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- (R 2 - M 2 ) 2 - (M 2 -N 2 ) 2 ) + 2ANR (4) 

- - (AM 3 + A(R 2 - M 2 )M 

2MNR y V ; 

- 4(Af 2 - N 2 )M), 

and cyclically in M, N, i?, pm,Pn, Pr- The dy- 
namical system (||) and (||) presents one critical 
point E with coordinates 



Pm = Pn=Pr = 0, M = N = R = M , (5) 

where Mo = \J~-^x- The energy of the critical point 
E is defined by 



2Er;rit A / . . • 



(G) 



Much of our understanding of nonlinear systems 
derives from linearization about critical points and 
the determination of existing invariant submani- 
folds, which are structures that actually organize 
the dynamics in phase space. The system presents 
a two dimensional invariant manifold of the dy- 
namics, defined by 



Pm =Pn = PR, M = N = R. 



(7) 



This invariant plane is actually the intersection 
of the two four dimensional invariant submani- 
folds defined by (M — N,Pm — Pn) and (iV = 
R,Pn — Pr)- The saddle-center-center critical 
point E belongs to the invariant plane, and from it 
emerge the separatrices S. The FRW singularity 
(M = N = R = 0,Pm — Pn = Pr = 0) is a degen- 
erate critical point of the dynamics on the invariant 
plane. The phase picture of the motion in the in- 
variant plane is given in Fig. 1, already expressed 
in canonical coordinates (x,p x ), to be introduced 




FIG. 1: Phase portrait of the invariant plane M = N = 
R,Pm = Pn = Pr, in the coordinates (cc,^) defined 
by (21) — (22). The separatrices S are characterized by 
the energy Eo = E cr it = • E is the critical point of 
saddle-center-center type, that belongs to the invariant 
plane (cf. text), x — 0,p x = is the FRW singularity. 



in Eqs. fl2l|)-([22|), that are variables defined on the 
invariant plane. Also, a straightforward analysis 
of the infinity of the phase space under consider- 
ation shows that it has two critical points in this 
region, corresponding to the DeSitter solution, one 
acting as an attractor (stable DeSitter configura- 
tion) and the other as a repeller (unstable deSitter 
configuration) for the dynamics at infinity. The 
scale factors M, N and R approach the DeSit- 
ter attractor as M = N = R ~ exp[t -\/A/3] and 
Pm = Pn — Pr. ~ exp[2ty^A/S\. It is easy to see 
that the DeSitter asymptotic configurations also 
belong to the invariant plane and that two of the 
separatrices S approach them for times going to 

±00. 

To proceed with the study of the phase space 
dynamics, let us now linearize the dynamical equa- 
tions (H), (Hp about the critical point E. We define 



X 
Y 
Z 

We obtain 



(M-M ),W = (p M -0), 
(N-M ),K = (p N -0), (8) 
(R-R o ),L = ( PR -0). 
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In what follows, without loss of generality, we fix 
A such that Mq = 1. The nature of the critical 
point E is determined by the characteristic poly- 
nomial associated with the linearization matrix 5". 
We obtain 



P(A)^(A 2 - i )(A 2 + 2) 2 , (9) 



These variables are conjugated to the pairs accord- 
ing to [X 1 , W] = 1, [K\ Y'] = 1, [L', Z'\ = 1, other 
Poisson brackets (PB) zero. The Hamiltonian ( |l2|) 
is separable, and we can immediately identify the 
following constants for the linearized motion 

E rotl = 3 -L'*+ l -Z>\ (14) 

Erot 2 = l -K' 2 +Y'\ (15) 

Q 1 = iY>Z> + ~K'L>, (16) 

Q 2 = L'Y' - Iz'K', (17) 



with roots 



A = ±i, A = ±iV2, 



(10) 



where the second pair has multiplicity two. The 
pair of real eigenvalues generates a saddle struc- 
ture, while the degenerate pair of imaginary eigen- 
values generate a double center structure. The 
analysis of the center structure will reveal a mani- 
fold of linearized unstable periodic orbits with the 
topology of a 3-sphere, as we will see. 

To display the structure of the linearized motion, 
we start by diagonalizing S with the use of a trans- 
formation matrix 3? whose columns are composed 
of six independent eigenvectors of S |Tq| , A judi- 
cious choice of 3? yields primed variables defined 
by the transformation 



X' 

y 
z 1 
w 

A" 
L' 



1 -{X + Y + Z), 

(X-Y), 

(X + Y -2Z), 

(W + K + L), 

\{W-K), 

1 
6 



(11) 



(W + K-2L). 



In these new variables, the quadratic Hamiltonian 
about E is expressed in the form 



H = {E crit - E ) + 1 (hv 12 - 6X' 2 



1 



yr. 



* L >2 \ z n 
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in the sense that they all have zero PB with the 
Hamiltonian (|l^). The first three constants ap- 
pear already as separable pieces in the Hamiltonian 
(|l2|). Ehyp corresponds to the energy associated 
with motion in the saddle sector, while E ro t 1 and 
E rot2 are the rotational energies associated with 
the motion in the linear part of the center-center 
manifold, namely, on stable periodic orbits in 2- 
dim tori. The remaining two constants are addi- 
tional symmetries that arise as consequence of the 
multiplicity two of the imaginary eigenvalues, and 
are associated to the fact that the linearized dy- 
namics in the center-center sector is that of a two- 
dimensional isotropic harmonic oscillator. They 
are not all independent but are related by 



4:E rotl E rot2 — 12Q 2 + 6Q 



(18) 



We are now ready to describe the topology of the 
general dynamics in a linear neighborhood of the 
critical point E. To describe all possible motions, 
the following situations must be taken into account 
in the linearized Hamiltonian (|l2|) . If Eh yp — two 
possibilities arise. First, we have W = X' = 0, 
implying that the motions are periodic orbits of 
the two-dimensional isotropic harmonic oscillator 



H=\-K>- 



Y' 



-L' 2 + -Z 1 ' 
2 3 

= {Ecrit — Eg 



(19) 



By a proper canonical reescaling of the variables in 
( [ill ) an d (19) it is easy to see that these constant 
energy surfaces are hyperspheres and that the con- 



(12) stants of motion Qi, Q 2 and Q3 — (E rotl — E, t 



rot 2 1 



satisfy the algebra of the three dimensional rota- 
tion group under the Poisson bracket operation, 
namely, 
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[QuQj] = e ijk Q k (20) 

The constant of motion Q\ considered as a gen- 
erator of infinitesimal contact transformations has 
a peculiar significance in characterizing the topol- 
ogy of the underlying group of the algebra (^0|) |17| . 
While Q2 generates infinitesimal rotations of the 
orbits, Qi generates infinitesimal changes in eccen- 
tricity. The action of Qi is to take an orbit - let 
us say nearly circular - and transforms it into an 
orbit of higher and higher eccentricity until it col- 
lapses into a straight line. Continued application 
of Q\ produces again an elliptic orbit, but now 
traversed in the opposite sense , so that it takes 
a 720° rotation to bring the orbit back into itself. 
The two-valuedness of the mapping arises from the 
fact that the orbits are oriented. Therefore we 
are led to the conclusion that the group generated 
by these constants of motion is the unitary uni- 
modular group and not the rotation group. It is 
therefore compatible to consider that the center- 
center manifold has indeed the topology of a three 
sphere S 3 . Now the separate conservation of E rotl 
and E rot . 2 (cf.Eq.(]I^)) allows us to show that the 
center-center manifold in the linear neighborhood 
of E is foliated by Clifford two-dimensional sur- 
faces in S 3 ]lq| , namely, two-tori ^se contained in 
the energy surface Eq = const. Such surfaces, as 
well as the S 3 manifold that contains them depend 
continuously on the parameter Eq- We remark 
that these two tori will have limiting configurations 
which are periodic orbits, whenever E rotl — or 
E ro t 2 = 0, and correspond to the case of maxi- 
mum eccentricity (for instance, a straight line in 
the plane {Y',Z')). 

The second possibility will be W = ±6X' , that 
defines the linear stable Vs and unstable Vjj mani- 
folds of the saddle-sector. Vs and Vu limit regions 
I (Ehyp < 0) and regions II (Eh yp > 0) of motion 
on hyperbolae which are solutions in the separable 
saddle-sector (cf. (Q) E hyp = \ (\W 2 - 6X' 2 ), 
(cf. (p^|)). Note that the saddle-sector depicts the 
structure of Fig. 1 in the neighborhood of E, with 
Vu and Vs tangent to the separatrices at E J|]. 
The direct product of $j_e with Vs and Vjj gener- 
ates, in the linear neighborhood of E, the structure 
of stable (^e x Vs) and unstable C^e x Vu) 3-dim 
tubes, which coalesce into the two-dimensional tori 
^Seq f° r times going to +00 or —00, respectively. 
The energy of any orbit on these tubes is the same 
as that of the orbits on the two tori 3b . These 
3-dim tubes are confined inside the 4-dim tubes 
which are the product of Vs and Vu by the S 3 
center-center manifold. It is obvious that the exis- 
tence of the S 3 center-center manifold is restricted 
to energy surfaces such that (E crit — Eq) > 0. The 
intersection of the center-center manifold with the 



motion 



away from 




the FRW 



FIG. 2: Stable and unstable 4-dimensional tubes em- 
anating from the center-center manifold S% . They 
are the nonlinear extension of S% (Qe ) x Vs and 
S% (Qe ) x Vu in the neighborhood of E. 

energy surface H = Eq is in general a three sphere 
parametrized with Eq, from which two pairs (one 
stable and one unstable) of 4-dim tubes emanate 
(cf. Fig. 2 and caption). More particularly, in the 
linear neighborhood of E (namely, for (E C rit — Eq) 
small), two pairs of 3-dim tubes emanate from the 
Clifford two-dimensional surface ^e contained in 
S 3 . For the case Eh yp 0, the motion is restricted 
on infinite tubes resulting from the direct product 
of the hyperbolae, lying in the regions I and II 
of Fig. 1, with Qe - A general orbit which visits 
the neighborhood of E belongs to the general case 
E h yp j= 0, E rotl ^ and E rot2 ^ 0. In this re- 
gion the orbits have an oscillatory approach to the 
3-dim linear tubes, the closer as Eh y p — > 0. 

We remark that the surface of the tubes consti- 
tutes a boundary for the general flow and, in the 
neighborhood of E, is defined by Eh yp — 0. De- 
pending on the sign of Eh yp the motion will be 
confined inside the 4-dim tube (for Eh yp < 0) and 
will correspond to a flow of motion distinct from 
the one confined outside the tube (for Eh yp > 0). 
This can be seen from Eq. (|l^) , that can be rewrit- 
ten 

E r oti + E rot2 — Eh yp = (E cr it — Eq). 

Let us consider, for instance, a set of initial con- 
ditions corresponding to initially expanding uni- 
verses. Two distinct flows will be associated with 
these initial conditions, depending whether they 
are contained inside the tube or outside the tube. 
A careful examination of Fig. 2 shows us that the 
flow corresponding to initial conditions inside the 
stable tube will reach a neighborhood of E and 
will return towards the FRW singularity inside the 
unstable tube (with E} lyp < 0), while the flow of 
orbits associated to initial conditions that are out- 
side the tube will reach the neighborhood of E and 
escapes towards the DcSittcr attractor along the 
exterior of the unstable tube that is directioned to 
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the right. Also the flow of the stable tube on the 
bottom right of Fig. 2 will escape towards the De- 
Sitter attractor along the interior of the unstable 
tube of the top right. We remark that the direction 
of the separatrices and hyperbolae in the linearized 
saddle-sector are the topological guide for the flow. 
Our main interest will rely on the stable and un- 
stable tubes that emerge from the neighborhood 
of the critical point E towards the FRW singular- 
ity and their homoclinic transversal crossings, that 
constitutes an invariant signature of chaos in the 
models. 

By continuity, the nonlinear extension of the 
center-center manifold will maintain the topology 
of the three-sphere S* 3 , but will not be decompos- 
able into E rotl and E rot2 so that now only the four 
dimensional tubes with topology Rx S 3 are mean- 
ingful for the nonlinear dynamics. 

The extension of the structure of the four di- 
mensional tubes away from the neighborhood of E 
are now to be examined, and our basic interest will 
reside in on the stable and unstable pair, S% x Vs 
and Sg o x Vjj, that leave the neighborhood of E 
and proceed towards the Friedmann-Robertson- 
Walker (FRW) singularity of the model. To this 
end let us consider the two dimensional invariant 
manifold of the dynamics, defined by Eqs.([5]). The 
two dimensional invariant plane is contained in a 
six dimensional phase space and it is obvious that, 
contrary to examples in lower dimensional systems, 
it does not separate the phase space in two disjoint 
parts. In fact the general motion about a neighbor- 
hood of the invariant plane is an oscillatory flow 
confined along the interior or the exterior of the 
four dimensional tubes R x S 3 , where the invari- 
ant plane (or more properly, one of the curves of 
the invariant plane) may be thought as a structure 
in the center of the tube, as we proceed to show. 
To see this, let us introduce the canonical coordi- 
nate transformation with the generating function 



1 

Pm = - 



NR 



3 {MNR)i 
1 NR 



P., 



1 N 

N Pv + & Pz 



AI 



M 



PN = 7T 



-Px 



3 (MNR)i 
1 NR 

PR 



N 2Py + R 2P z 



(22) 



2MN 



3 (MNR) 3 * 



R? 



-Pz 



It is worth remarking that the linearization of this 
canonical transformation about the critical point 
E yields exactly the linear transformation (|TT|) , 
and that the variables (y,p y , z,p z ) correspond to 
the primed variables (K 1 , Y', L', Z') defined on the 
S 3 center-center manifold about a linear neighbor- 
hood of the critical point. The variable x is ob- 
viously the average scale factor of the model. In 
these new variables, the full Hamiltonian (|2|) as- 
sumes the form 



H 



x 
yzi 



" 24x Px 
x 1 

~2^~ 2 
yx 



]f_ 2 _ 3f! 2 
2x 3Py 2x 3Pz 



xz$y 



1 



rXZ* 



XZ^ 



2 3 



2 Ax 3 = 2E , (23) 



and the equations of the invariant plane reduce to 



V =1,Z=1, Py = 0=Pz- 



(24) 



It is clear that {x,p x ) are variables defined on 
the invariant plane, and now the expansion of the 
Hamiltonian ( p3| ) about the neighborhood of the 
invariant plane can be easily implemented, pro- 
ducing a linearized Hamiltonian parametrized by 
the variables (x(t),p x (t)) describing the curves in 
the invariant plane, in analogy with the way we ex- 
pand a dynamical system about a periodic orbit. 
We obtain 



i M MN H = -^—pl. + — - 2Ax 3 ^rvl 

G=(MNR)3 Px + ^p y + — Pz 2Ax lx 2 2x 3 ^ 2x 3 ^ 



-x(y - l) 2 - \x{z - l) 2 = 2£ (25) 



with dynamical equations 

where p x , p y and p z are the new momenta, result- 
ing in . i 

(Sp y ) = 2x5y 

x =(MNR)i,y=^,z = ^- (21) {5z) = ~5p z (26) 



and 
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where Sy = y — 1, <5z = z — 1, (5p y = p y — and 
<5pz = p z — 0. The linearization matrix in ( |26| ) has 
imaginary eigenvalues, both with multiplicity two, 
given by A = ±^5, so that in the neighborhood 

X{t) 

of the invariant plane we have only elliptic modes, 
namely, the motion is oscillatory about the invari- 
ant plane with increasing frequency as the orbits 
approach x = 0. Close to the critical point E, the 
motion in the invariant plane is slow, i.e., (x,p' x ) is 
small. Following an orbit in this plane with dimin- 
ishing x, its phase space speed increases but not 
faster than the frequency of the surrounding mo- 
tion in the other four variables. We are thus justi- 
fied in deriving a qualitative picture for the motion 
near the invariant plane from an adiabatic approx- 
imation for the relevant orbit manifolds. That is, 
we will assume that the overall geometry of the 
four dimensional stable and unstable tubes that 
are asymptotic to the S 3 center-center manifold 
can be approximated by the adiabatic deformation 
of this sphere as x decreases away from the criti- 
cal point. In this way, we have in the six dimen- 
sional phase space the four dimensional tubes of 
motion R x S 3 with the two dimensional invariant 
plane as the structure at its center, about which 
are the oscillatory degrees of freedom (y,p y , z,p z ). 
We will use this fact to show that the motion in the 
invariant plane guides the four dimensional tubes 
of motion inducing inevitably the crossing of the 
unstable tube with the stable one in the neighbor- 
hood of (a; = p x — 0). Also the eigenvectors of 
the matrix in (Eq) can give us an idea of the be- 
haviour of the flow in the tubes as x diminishes, 
although in the neighborhood of x — the linear 
expansion (|25|)-(|26|) is no longer valid. Associated 
with the eigenvalue A = we may choose the 
two independent eigenvectors 



ei 



C2 



1 2x 



V2 



0,0 



x i 



1 2x 
0,0,-, 



x 3iV2. 
and for the eigenvalue A = — 



e 4 



iV2 
~2x~' 



£,0,0 



0,0, — ,x 

2x 



From the above form of the eigenvectors, we can 
draw a series of important informations about the 
behaviour of the four dimensional tubes as they 



approach the singularity x = 0. As a: decreases, 
the sections x = const, of the four dimensional 
tubes about the invariant plane stretches in two 
directions and contracts in other two, while main- 
taining their symplectic invariants constant. Also 
as x decreases the oscillations have increasing in- 
stantaneous frequency cx This appears as a 
characteristic of the dynamical approach to the 
singularity even for the extension of the tubes to 
a nonlinear neighborhood of the invariant plane. 
There is a further fundamental fact that results 
from the ^-dependence of the above eigenvectors: 
let us consider for instance the sections x = 1 of 
the tubes which emanate from the two-tori $se 
defined in the linear neighborhood of the critical 
point E. Their image projected on the planes 
(y,p y ), (z,Pz) can be rescaled to small circles that 
deforms into ellipses, as x is infinitesimally dimin- 
ished; the ellipse is rotated infinitesimally coun- 
terclockwise if the section corresponds to the in- 
ferior branch of the curves in the invariant plane 
(p x < 0) or clockwise if it corresponds to the su- 
perior branch (p x > 0). This will be fundamental 
to garantee the transversal crossing of the tubes in 
the neighborhood of x — 0, as we will discuss. We 
must however remark that this conclusion is taken 
in the adiabatic approximation, namely, consider- 
ing that the variation of x with time is adiabatic; 
this certainly does not occur near the singularity 
but wc assume that no drastic change in the dy- 
namics will alter this behaviour for the actual dy- 
namical tubes. We illustrate this numerically in 
Fig. 3 by depicting the sections x — 0.5 projected 
on the planes (y,p y ). Figs. 4 show a numerical 
illustration of the unstable and stable 4-dim tubes 
emanating from the neighborhood of E towards 
the FRW singularity. We note that, for E cr i t — Eq 
small, the motion proceeds towards the singular- 
ity along a tube whose projection on the invari- 
ant plane (x,p x ) shadows the invariant curve cor- 
responding to H = Eq, along the upper branch 
and the lower branch (cf.Fig. 4(b)). 

We proceed now to show the homoclinic crossing 
of the stable and unstable tubes at the neighbor- 
hood of the singularity (x = 0,p x = 0), a phe- 
nomenon analogous to Poincare's homoclinic tan- 
gle, and source of chaos in the model[^0|. The 
four following points are essential for this. First, 
in the five dimensional energy surface H = Eq, 
the four dimensional tubes of motion R x S 3 ac- 
tually divide the space, separating the dynamics 
outside from the dynamics inside the tubes. Sec- 
ond, in the canonical variables {x,p Xl y,p y ,z,p z ) 
introduced above, we have seen that (x,p x ) are the 
variables on the invariant plane, whose phase space 
picture is depicted in Fig. 1, and the remaining de- 
grees of freedom (y,p y , z,p z ) correspond to elliptic 
modes of motion (cf. analysis in the neighborhood 
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of the invariant plane). Actually, the tubes have 
the two dimensional invariant plane as the struc- 
ture at its center (more properly, one of the curves 
of the invariant plane corresponding to H = Eq), 
about which the flow with the oscillatory degrees 
of freedom (y,p y , z,p z ) proceeds and in fact the 
motion in the invariant plane guides the flow, in- 
ducing inevitably the crossing of one of the four- 
dimensional unstable tube with the stable one in 
a neighborhood of (x = 0,p x = 0). Third, we 
note that the invariant plane cannot leave the in- 
terior of the four dimensional tubes since they are 
separation surfaces in the five dimensional mani- 
fold H — Eo, and fourth the property of clock- 
wise rotation or counterclockwise rotation of the 
sections of the stable and unstable tubes, respec- 
tively, in their approach to the FRW singularity. 
Summing up, in the five dimensional energy sur- 
face H = E , the model presents the structure of 
four dimensional tubes of orbits, emanating from a 
three sphere along the two saddle directions about 
the critical point saddle-center-center E, one sta- 
ble and one unstable (the stable one correspond- 
ing to initially expanding universes). These tubes 
have the invariant plane curve as the structure at 
its center and the motion in the invariant plane 
guides the oscillatory flow towards the singular- 
ity (a; = 0,p x =0), inducing inevitably the cross- 
ing of the unstable tube with the stable one in a 
neighborhood of the singularity; that the crossing 
is not just a point comes from the conservation 
of sympletic areas in Hamiltonian dynamics, and 
the transversality of the crossing comes from the 
clockwise or counterclockwise rotation of the sec- 
tions (y,p y ) and (z,p z ) of the stable and unstable 
tubes, respectively. If we consider the transver- 
sal crossing in a section, say x — const., it is not 
difficult to see that the intersection is a S 2 mani- 
fold. Therefore the intersection manifold will be a 
3-dim tube of flow (with topology R x S 2 ) which 
is contained both in the J^-dim stable tube and in 
the 4- dim unstable tube, and homoclinic to the S 3 
center-center manifold. It is the equivalent of a 
homoclinic 1-dim orbit in lower dimensional cases. 

Typically, if the 4-dim tubes intersect transver- 
sally once they will intersect each other an infinite 
number of times producing an infinite set of homo- 
clinic orbits which are actually topological 3-dim 
tubes R x S 2 . The homoclinic 3-dim tube, which 
is bi-asymptotic to the S 3 manifold provides the 
mechanism for stretching and contraction, giving 
origin to the homoclinic tangle, which is an in- 
variant signature of chaos in the model JnJ. We 
just mention that the dynamics near homoclinic 
orbits is very complex, forming well-known horse- 
shoe structures (cf. || , @ , Q and [||, and the 
bibliography therein) . 

Let us consider the first intersection of the tubes; 




0.9999 1 1.0001 

V 



FIG. 3: Sections of the stable and unstable tubes pro- 
jected on the plane (y,p y ). The dotted circle is the 
projection of the 2-dim torus Ss_e defined in the lin- 
ear neighborhood of E. The ellipses are the sections 
for x = 0.5 of the tubes emanating from ^5e - Note 
the distortion of the circle into the ellipses and the 
clockwise/counterclockwise rotation of the sections of 
the unstable/stable 3-dim tubes as x diminishes (from 
x — 1 to x = 0.5) towards the FRW singularity. 



a part of the orbits inside the unstable tube will 
enter in the interior of the stable tube and the 
flow will proceed along the stable tube towards the 
neighborhood of the critical point E, from where 
it will re-enter the unstable tube and proceeds to- 
wards the FRW singularity and by a new intersec- 
tion a part of these orbits will again enter the stable 
tube and proceed back towards the neighborhood 
of the critical point E, and so on, producing an in- 
finite recurrence of the motion for a class of orbits, 
which will be periodic orbits of very long periods 
and bounded oscillatory orbits. Another part of 
the orbits inside the unstable tube will flow along 
the exterior of the stable tube towards the neigh- 
borhood of the critical point E, and afterwards 
will escape towards the DeSitter attractor at infin- 
ity along the exterior of the unstable tube of the 
second pair (cf. Fig. 2). This for one single in- 
tersection. The same pattern will be reproduced 
for each of the infinite intersections of the unsta- 
ble tube with the stable one. In this process, the 
surface of the tubes, which is in fact a boundary 
for distinct types of flows, will become more and 
more stretched and folded, resulting in an intricate 
structure which will have some important physical 
consequences for the long time behaviour of the 
orbits. Indeed, basins of initial conditions for ini- 
tially expanding universes have a fractal boundary 
associated with recollapse or escape into the De- 
Sitter configuration, as we show soon. 

We must comment that the chaotic aspect of the 
dynamics described here depends crucially on the 
presence of the cosmological constant that engen- 
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Fig. 4(a 




Fig. 4(b) 



FIG. 4: (a) Numerical illustration of the stable and 
unstable tubes emanating from the neighborhood of 
the critical point E towards the FRW singularity, pro- 
jected on the submanifold (x,y,p x ), for E cri t — Eo = 
10 -8 . (b) Projection of the same stable and unstable 
tubes of (a) on the invariant plane {x,p x ). Note that 
the projection "shadows" the orbits on the invariant 
plane corresponding to H — Eo- The insets show the 
sections x = 0.2 projected on the (y,p y ) plane. 



ders the structure of a saddle-center-center critical 
point in phase space together with the homoclinic 
four dimensional tubes, the transversal crossings 
of which (in a neighborhood of the FRW singular- 
ity) characterize chaos definitely and unambigu- 
ously in the model. If we restrict ourselves to 
the energy surface Eq — we have the Mixmas- 
ter universe with a cosmological constant which 
presents the same structure of chaos. For the lim- 
iting case (Eo = 0,A = 0), the absence of the 
critical point appears to eliminate the set of ho- 
moclinic orbits and periodic orbits of arbitrarily 
large periods, whose recurrence engenders chaos in 
the models. In this instance, our approach does 
not allow to invariantly characterize chaos in the 



Mixmaster model. 

Our system is open (noncompact) with two def- 
inite asymptotic exits, namely, collapse and es- 
cape to the DeSitter configuration. The code es- 
cape/collapse defines basin boundaries in the ini- 
tial conditions set; these boundaries are associated 
precisely to the surface of the homoclinic tubes in 
the set, and correspond to the initial conditions 
for the homoclinic intersection manifold contained 
simultaneously in the stable and in the unstable 4- 
dim tubes. Together with the countable set of pe- 
riodic orbits of arbitrarily large periods that exist 
in the neighborhood of each homoclinic orbit and 
that have the homoclinic orbit as an accumulation 
set [ 2jJ p3] , they constitute the set of orbits which 
neither escape or collapse. This set of bounded 
orbits that are in the boundary of collapse and es- 
cape constitutes , in this way, the strange repellor. 
The fractal dimension of the basin boundary sets 
are calculated below by a box-counting method. 



FIG. 5: Plot of the scaling law / cx e a in log- 
log scale, where e is the uncertainty radius about 
5,000 initial conditions taken inside the set 23 con- 
structed about the point (x = 0.4, y = l, z = l,p x = 
2.036467455281656, p a = 0,p z = 0) on the separatrix 
of the invariant plane (x,p x ). f is the fraction of 23 of 
uncertain initial conditions with respect to the uncer- 
tainty code collapse/escape. The set 73 is a hypercube 
of edges 10 -4 about the mentioned point, in the energy 
surface E = 0.9999999683772234. The best fit of the 
small linear region renders (a) a ~ 0.321. 

In other words, our system has the characteristic 
of a chaotic scattering system, having the saddle- 
center-center as a chaotic scatterer, with two abso- 
lute outcomes consisting of (i) escape into inflation 
(the DeSitter attractor) or (ii) recollapse to the sin- 
gularity, for initially expanding universes. The De- 
Sitter attractor defines a way out from the initial 
singularity to the inflationary phase but this exit 
to inflation is completely chaotic, namely, small 
fluctuations in initial conditions may cause the 
universe to change its asymptotic behaviour from 
recollapse to escape into inflation and vice-versa. 
This is due to the presence of chaotic sets in the 
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phase space of the system, as a consequence of the 
structure of homoclinic tubes and their transversal 
crossings, which provides an invariant characteri- 
zation of chaos in the model. One of the ques- 
tions to be examined here is the characterization 
of sets of initial conditions for which the DeSit- 
ter attractor is attained and establish the char- 
acter of a chaotic scattering system. To see this 
chaotic exit to inflation, let us consider initial con- 
ditions sets for initially expanding universes in the 
energy surface E such that the orbits visit a lin- 
ear neighborhood of the saddle-center-center crit- 
ical point E. In this neighborhood, the Hamilto- 
nian is separable according to (12), namely, H = 
(Ecrit - E a ) + E hyp - E rotl - E rot2 = 0. The par- 
tition of the energy (E cr n — Eq) into the energies 
E^p and E rotl + E rot2 of motion about the criti- 
cal point will determine the outcome of the oscil- 
latory regime about E into collapse or escape to 
inflation (DeSitter attractor) whether Eh yp < or 
Ehyp > 0, respectively. However the nonintegrabil- 
ity of the system, with the consequent homoclinic 
crossing of the tubes, will cause this partition to 
be chaotic in general, namely, given an arbitrary 
initial condition of energy Eq we are no longer 
able to foretell in which of regions I or II about 
the saddle-center the orbit will land when it ap- 
proaches E. Since Eh yp = is a limiting case, for 
initially expanding universes, the motion inside the 
tubes corresponds to orbits that will recollapse af- 
ter reaching the neighborhood of E, while the mo- 
tion outside the tubes corresponds to orbits that 
will escape into the inflationary phase (towards the 
DeSitter attractor at infinity). In other words, the 
intrincate crossing and merging of the tubes pro- 
duces chaotic sets in the phase space of the model, 
in particular establishing fractal basin boundaries 
associated with the chaotic exit to inflation. To 



illustrate this, we finally present a numerical ex- 
periment where we obtain a measure of the frac- 
tal dimension of basin boundaries in sets of ini- 
tial conditions by using the box-counting method 
with the uncertainty code collapse/escape into in- 
flation. In this numerical experiment we evalu- 
ate the fractal dimension of portions V of phase 
space, about a point of the separatrix S on the 
invariant plane. The box-counting method, due 
to Ott and collaborators 23 , consists of deter- 



mining the uncertainty exponent a appearing in 
the scaling law / oc e", where e is the uncertainty 
radius about 5,000 initial conditions taken inside 
the sets T>, and / is the fraction of T> of uncer- 
tain initial conditions with the uncertainty code 
collapse/escape into inflation. The uncertainty ex- 
ponent is related to the fractal dimension d of the 
fractal basin boundary by d = N — a, where N = 5 
is the phase space energy surface dimensionality. 
In our numerical experiment we selected a set T> 
of homogeneously and randomly distributed initial 
conditions about the point [x = 0.4, y = l,z = 
l,p x = 2.036467455281656, p y = 0,p z = 0) on the 
separatrix S, such that the orbits visit a neighbor- 
hood of the saddle-center-center E before collapse 
or escape into inflation. In Fig. 5 we depict the 
scaling law / cx e a in which the small linear portion 
lies between the saturation region for large e and 
noise for very small e. The following value for the 
fractal dimension of the set V is found, d w 4.627, 
confirming the chaotic nature of these sets. 
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